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CN ■ Observations over the past two solar cycles show a highly irregular pattern 

■ of occurrence for major solar flares, 7-ray events, and solar energetic particle 
(SEP) fiuences. Such phenomena do not appear to follow the direct indices of 
solar magnetic activity, such as the sunspot number. I show that this results from 
non-Poisson occurrence for the most energetic events. This Letter also points out 

CLc a particularly striking example of this irregularity in a comparison between the 

O ■ declining phases of the recent two solar cycles (1993-1995 and 2004-2006, respec- 

. tively) and traces it through the radiated energies of the flares, the associated 

. SEP fiuences, and the sunspot areas. These factors suggest that processes in the 

solar interior involved with the supply of magnetic flux up to the surface of the 
^ . Sun have strong correlations in space and time, leading to a complex occurrence 

QQ . pattern that is presently unpredictable on time scales longer than active-region 

y—^ . lifetimes (weeks) and not correlated well with the solar cycle itself. 

^ . Subject headings: Flares - Solar Energetic Particles - X-rays - Sunspots 

i> : 
o ■ 

> ■ 

■ 1. Introduction 

. 5t , It has long been known anec dotally that highly en ergetic solar events do not strictly 



follow the solar sunspot cycle (e.g., iGarcia fc Dryerlll987l ). The fact that we only have a few 



cycles of modern data has made it difficult to describe this discrepancy quantitatively, espe 
cially in view of the small numbers of the most energetic events. The fossi l records typically 



do no t have enough time resolution for to overcome these problems (but see iMcCracken et al. 



2OOII ). These most energetic events include some of the most geoeffective ones, so there we 
have a clear practical reason for studying their occurrence patterns - we would like to predict 
the occurrence of a major event. 

The most energetic events also represent the extreme limit of the mechanism that stores 
energy in the solar corona. In the consensus view magnetic energy builds up gradually in the 
corona as a result of stresses imposed from below. The stressed coronal field then relaxes. 
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by unknown processes, to produce a flare and/or coronal mass ejection (CME). The energy 
appears to a rrive in the cor ona as the result of buoyant motions of current-carrying flux 
systems (e.g.. lSchrijverll2007l ) rather than by the twisting of the coronal field by photospheric 
surface flows as often assumed in numerical simulations. The patterns therefore reflect the 
persistence of t he flux-emergence process , which is known to display coherence in both space 



and time (e.g., iKnaack fc Stenflol l2005l ). and ultimat ely must be attr ibuted to the solar 



dynamo and other processes in the solar interior (e.g., RuzmaikinI 1 1 9 Qsl ) . 



Flare occurren ce apparently follows a nonstationary Poisson dis tribution with time- 
varying mean rates (lBieseckerlll994l : IWheatlandll2000l : iMoon et al.ll200ll ) and a clearly power- 
law dependence on event "size," where this conceptually reflects total event energy but in 
practice often refers to an observational parameter such as peak X-ray luminosity (e.g.,|Drake 
197ll : lHudsonlll99ll ) . Many studies have shown that flare occurrence follows a flat power-law 
relationship, d{logN) / d{logE) = —a, with a < 2. There are suggested weak dependences 



of the exponent on the phase in the solar cycle (jBailll993l : IWheatland fc Litvinenko 



by active region (IKucera et al.l 119971 ). and for from star to star (e.g., IShakhovskaya 



20021) 



19891). 



Such a flat distribution requires a high-energy cutoff to conserve energy, but there is no clear 
evidence for such a cutoff yet. 

The more energetic the flare, the more likely the occurrence of a CME, although in a 



de La Beauiardiere et al. 



lYashiro et al. 



few c ases an X-class flare will not have a CME association (e.g 
1995). For weaker flares, associated CMEs occur much less frequently (e.g 
2006,) . The CME distribution must therefore devia te from the flare power la w at low event 
energies, possibly not following a power law at all (I Jackson fc Howard! Il993l ). Interestingly, 
solar e nergetic particle fluences do follow a power law, but a significant ly flatter one th an the 
flares (Ivan HoUebeke et al.l 1 19751 : iGabriel fc Feynmaru Il996l ) ; see also iHudsoru (119781 ) . The 
occurrence of solar energetic particles (SEPs) might otherwise be expected to re flect the 
CME distribution, beca use CME-driven shocks are known to accelerate SEPs (e.g.. iReames 
I999I : Iciiver et al.ll2004h . 



In this Letter we report a large specific variation in X-class flare occurrence rate that we 
trace through similar patterns in SEP fluences and in sunspot areas. This juxtaposition is 
consistent with the interpretation of flare occurrence with Biesecker's variable-rate Poisson 
process, although the small numbers of the most energetic flares means that this interpre- 
tation is only weakly grounded in this context. We instead suggest an origin in correlations 
of solar interior magnetism on time scales longer than about one rotation period, whose 
existence will strongly limit flare prediction on these time scales until the interior dynamics 
is better understood. 
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2. X-class Flares 

An X-class flare corresponds to a peak flux of 10^^ W/m^ in the GOES standard 2- 
8A passband. Such events he at the upper end of the occurrence energy distribution function 
of all flares, and may differ in their temporal occurrence because of the requirement for 
an upper energy cutoff - because of this, one cannot assume that the energy distribution 
continues to have the same power-law form as the flaring rate changes. Their small numbers 
(about 125 in the past solar cycle, from 1996 through 2006) make statistical analyses difficult, 
and in fact the more energetic of these events may saturate the detectors, which tends to 
diminish the quality of the statistics. 

The declining phases of the past two solar cycles have shown a striking discrepancy in 
the occurrence of X-class flares. This got attent ion because of the RHESSl observations 



of 7-ray flares in 2003-2005 (e.g.. IShih et al.ll2006l ): such events typical ly correspond to the 



X-class flares, and RHESSl observed several remarkable examples (e.g.. lShare et al.ll2004l ) in 
its inaugural years 2002 and 2003. The expectation for the years 2004-2006, if based on the 
previous-cycle years of approximately 1993-1995, would have been zero further events - not 
a single X-class flare occurred during these three late years of the previous cycle, alth ough 



one old-cycle event did occur in 1996 (IKosovichev fc Zharkovalll998 



Hudson et al.lll998h . To 



our surprise as many as 34 X-class flares occurred over 2004-2006, though not all observable 
as 7-ray events from RHESSl because of its orbital eclipse cycle. See Figured] for the data, 
all of which were obtained from Web resources maintained by NOAa|^ 

Figure [1] shows three cycles of X-class flare occurrence, highlighting the discrepant be- 
havior in the decaying phases of Cycles 21, 22 and 23. The difference in occurrence of 
energetic events between the latter two epochs is highly signiflcant; for a guide to signifi- 
cance we can use a Poisson distribution based on the number of unique active regions in 
the years 2004-2006 (11 unique regions, for an average of about 3 X-class flares per region). 
Computing the Poisson probability of one event in the earlier epoch (the 1996 flare) relative 
to the number of unique regions of the later epoch, we find a likelihood of <0.02%. This 
conservatively confirms the obvious inference from the Figure, namely that the X-class event 
numbers are highly discrepant and that the occurrence of such major energetic events has 
shown much greater variation than the sunspot number itself. Cycle 21, on the other hand, 
showed an intermediate number of events (15 X-class fiares, from 9 unique regions) and does 
not appear discrepant. 



^ ht tp : / / www . ngdc . noaa. gov /stp/ S OL AR 
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3. Solar Energetic Particles 



The striking difference shown by the X-class ffare occurrence between the past two cycl e 
declining phases also shows up strongly in the SEP ffuences (Figure [2l from iReedyl l2006l ) . 
This would be expected because of the strong c orrelation between X -class flare occurrence 
and CME occurrence, as documented recently by lYashiro et al.l (120061 ) . The declining phases 
of the two recent cycles, comparing (for example) 1994 with 2005 in Figure [21 clearly differ 
significantly. 

The identification of flare activity with SEP fluxes might seem inconsistent with the 
theory of particle accele ration by CME-driven shocks, rather than flares per se (e.g. jReames 
19991 : ICliver et al.ll2004l ). and frequent assertions of the i ndependence o f CME and flare oc- 
currence. This becomes understandable from the work of lYashiro et al.l (120061 ) . who confirm 
the well-known strong association of CMEs with the most energetic flares. The discrepancy 
in the numbers of the most energetic events between the two recent cycle declining phases 
can thus be traced in flare, CME, and SEP occurrence patterns. We discuss the significance 
of this finding in Section [5] but first investigate whether or not this occurrence discrepancy 
can also be detected in sunspot area statistics. 



4. Sunspot areas 



The plot in Figure [3] shows data obtained from the tabulations of sunspot group area by 
the SOON^ stations. A large fraction of the tabulated data have been used, typically from 
three or more stations for each day, but with rejection of a small number of outliers and also 
the measurements with quality values below 3 (the range is 1-5; see the NOAA Web site for 
details of the SOON sunspot data). The solid line in the plot shows the mean of the maxima 
of the daily areas for individual groups, in millionths of the hemisphere (the customary unit). 
This shows a time variation significantly distinct from that of the number of groups (dotted 
line) which roughly tracks the sunspot number. The larger values of mean areas during the 
decay phase of Cycle 23 (2004-2006) shows that the distribution function of sunspot group 
areas favored larger spots than during the corresponding interval in Cycle 22 (1993-1995). 
This asymmetry coincides with the asymmetry noted above in X-class flare occurrence and 
in SEP production. 
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Discussion 



Major energetic solar events do not closely track the solar cycle as a source of the 
slow variation under the dominant Poisson statistics. Indeed, the "Bayesian blocks" of 

time scales for Poisson behavior obtained by other methods (e.g., 



Wheatland 


( 


2000) or 


Gallagher et al. 


2002) 



events (on the order of one event per month over 1996-2006). We conclude that other physics 
dictates the occurrence patterns of the most energetic events, for which at most a few may 
occur in a given active region. The underlying cause of the Poisson behavior for the less 
energetic events should be found in the physics of energy buildup and release in the corona. 
The occurrence of the most energetic events presumably has more to do with the broad-band 
co herence of solar magnetic activity on large scales in both space and time, as discussed 



by iKnaack fc Stenfid (120051 ) in terms of "intermittent oscillations" revealed by spherical- 
harmonic expansions of synoptic magnetogram data. Examples of broad-band correlations 
would include the butterfly diagram and the presence of "active longitudes" where active 
regions may occur repeatedly. We can also note the remarkable eruption of three distinct 
active regions in October 2003, each producing X-class flares, and with distinct active regions 
in both hemispheres. Such a sudden and widespread surge of activity is certainly remarkable, 
even though noted here only a posteriori. 



Magnetic flux emergence leads directly to flare activity (e.g. ISchrijven 120071 ). and the 
occurrence of multiple major flares in a given active region therefore points to a persistence 
in the pattern of flux emergence. This persistence seems to be required to explain the 
occurrence of homologous flares, since we believe that extracting the energy from stressed 
coronal magn etic fields require s thei r irreversible restructuring, for example by magnetic 



reconnection. 



Nitta fc HudsonI (120011 ) show that this persistence can result in homologous 
CMEs in association with impulsive X-class flares. For reasons currently unknown, the 
strongest flux emergence, leading to the most energetic solar events, does not follow the 
relatively smooth pattern of flux emergence that defines the solar cycle and the occurrence 
patterns of less-energetic events. 

The striking variability in the occurrence of energetic events described in this paper 
might correspond to a modulation of the event rate near the upper limit on flare energy. 
Such a cutoff is required by the non-convergence of the flat occurrence power law of solar 
flares. The existence of a cutoff in particle fluences is already well-established from the 
fossil records, which have the advantage of extending over longer periods of time and thus 
of capturing the rarer extremely energe tic events. The ^'^C record suggests a maximum 



SEP fluence of some 10 ^'^ protons cm ^ (ILingenfelter &: Hudsonlll980l) a n d fossil cosmic-ray 



records over longer time scales agree well with this limit (IReedyl Il996l ). iMcCracken et al. 
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( 120011 ) set the cutoff at about 6 x 10^ protons cm~^ (omnidirectional ffuence) at >30 MeV 
based upon nitrate concentrations in Greenland ice cores. This proxy has the advantage that 
it overlaps the historical record. 

The SEP cutoff ffuence co rresponds roughly to the largest X-ray ffare ffuxes, of class XIO 
( iLingenfelter fc Hudson! Il980l ). Observing an analogous cutoff in the X-ray ffuxes (or other 
measures of ffare energy) is difficult, however, both because of the rarity of the most energetic 
events and also because they tend to cause detector problems that make it difficult to obtain 
precise photometry (the GOE^I photometers themselves saturate at about this level). Such 
a cutoff i n X-ray ffare sta tistics, which best reffect total ffare energy, has not yet been 
reported. iNita et al.l (120021 ) actually do observe an upper cutoff in radio burst magnitudes, 
in a comprehensive study, but they also note calibration difficulties and other factors t hat 
may contribute to this. The SEP fluxes have a "streaming ffux limit" (e.g.. lReameslll999l ). so 



the agreement of the SEP cutoff with the presently-observed maximum in the GOES event 
energies may be fortuitous. 

Does any index of ffare magnitude show a similar high-energy limit? The soft X-ray 
photometry from GOES provides the most stable long-term database of ffare magnitudes, 
and we have analyzed it to answer this question. Figure H] shows the distribution of M- and 
X-class ffares for the period from September, 1975, through January, 2007. This consists of 
5,637 M events, 424 X events, and 22 "super" events above XIO (numbers inclusive of Ml.O, 
XI. 0, and XIO.O). We do not show the super events in the Figure b e cause of distortion due 
to saturation. The maximum- likelihood method of ICrawford et al.l (Il970l ). independent of 
binning, gives a fit over the M-X range of dn/dS = 5520 x 5'-2-i93±o.oi5 gyents per unit 
X-class interval, the differential distribution. This distribution predicts 24.6 super-events, 
whereas 22 were actually observed. Within errors, there is thus no downward break. The 
fit over the M-X range given here is slightly steeper than expected, probably because of the 
lack of background subtraction in the reported event magnitu des. The ffare e nergy upper 
limit must therefore be significantly above XIO - as noted by ISchaefer et al.l (120001 ) . solar 
super-events, were any to have occurred, ought to have been detected by solar astronomers 
within the historical era. 

Resolving this question - at what point does the ffare energy distribution steepen? - 
would provide a important c lue for studen t s of th e generation of solar magnetic ffux and its 
delivery to the photosphere. iKucera et al.l (119971 ) interestingly suggest that a cutoff may be 
observable directly in event distributions for smaller active regions, at lower event energies. 
Thus the hypothetical cutoff in X-ray ffare magnitudes might reffect the downturn in active- 



^ Geostationary Operational Environmental Satellite 
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region areas expected from the log-normal distribution noted for sunspot areas (iBogdan et al. 



19881 ). The result regarding mean areas (Figure [3]) conflicts with the stability of the spot area 



distribution noted by Bogdan et al., but this may reflect the differing time scales studied . 
The existence of the needed cutoff in the distribution has been anticipated by lMuUanl (Il975l ). 
who suggested relating the maximum energy of a stellar flare with the scale lengths present 
in the convection zone of the star. 



6. Conclusions 

We have shown, based on the decay phases of Solar Cycles 22 and 23, an unexpected 
example of large-amplitude variations in the occurrence of the most energetic solar events. 
We could also trace this pattern in SEP fluxes and in sunspot group areas. These most 
energetic events (GOES XI or greater) do not follow the usual Poisson statistics with mean 
rates that govern lesser flares with shorter waiting times. The waiting times for the most 
energetic events indeed often exceed the active-region lifetimes, or the solar rotation period. 
Their statistics therefore reflect physics unrelated to coronal energy buildup and the mean 
flaring rate for a given active region. We suggest that solar interior dynamics dictates the 
pattern of occurrence of the most energetic events, rather than the coronal development. 

This draniatic variability reduces the predictability of major hazards in space (e.g.. 
Smith fc Scald 120071 ). since it is clear that a variable-rate Poisson distribution following the 
solar cycle as defined by a smooth sunspot number will not suffice. W orse yet, the flatness 
of the particle fluence distribu tion - which has an index of 1.2-1.4 (Ivan HoUebeke et al. 



1975 



Gabriel &: Fevnmaru 119961 ). flatter still than the flare energy distribution at about 1.8 
(e.g.. iHudsonI 1991 ) - means that individual events will dominate the total X-ray and 7-ray 
fluences. At present such events are basically unpredictable on time scales longer than a few 
days. 
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1975 1980 1985 1990 1995 2000 2005 

Year 

Fig. 1. — X-class flare numbers by year from 1 September 1975 tlirougfi 31 Jamiary 2007. 
Points sliown as diamonds are tlie years of the solar-cycle declining phases, defined here as 
1983-1985, 1993-1995, and 2004-2006. The corresponding numbers of X-class flares are 15, 
zero, and 34 respectively. 
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Fig. 2. — Solar energetic particle (SEP) event occurrences for 1954-present (excluding the 
events of 2006 December), show n as dashed vertical lines for >10 MeV threshold and solid 
vertical lines for >30 MeV (from lReedyll2006l ). The background curve is the sunspot number 
in monthly bins. Note the large fluences around 2005, and the negligible fiuences one cycle 
earlier around 1994. 
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Fig. 3. — Sunspot and flare behavior during Cycles 22 and 23. Dotted line, the annual 
numbers of sunspot groups; solid line, 2 x the mean peak areas of the groups (see text). 
Histogram, the numbers of X-class flares x 8. The vertical dashed lines mark the two 
declining-phase epochs studied in this paper. Data from the SOON network via NOAA. 
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Fig. 4. — Distribution of GOES 1-8A peak fluxes for the interval September, 1975, through 
January, 2007, for the M and X-class events (discarding the 22 "super-flare" occurrences 
above XIO). The d ashed lines shows a fit using the maximum likelihood method of 



Crawford et al.l (Il970l ). which does not require binning. The binning shown is 0.1 X units 
for the M flares, and 1 X unit for the X flares (where XI corresponds to 10~^ W/m^ peak 
soft X-ray flux). This fit predicts the observed number of super-flares within errors, giving 
a lower limit on the break energy. 



